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Abstract 

In this paper, I aim to study free oscillations of a system of oscilla- 
tors in more than one dimensions in the absence of damping. The basic 
approach lies in decoupling the motion in the individual perpendicular 
directions. Once the equations are decoupled, the existent techniques 
of Normal mode expansion for 1-dimensional oscillators are used to 
solve for the equations of motion. I also study the motion of a driven 
system of oscillators in higher dimensions in the presence of a velocity 
dependent damping force. 



1 Introduction 

Motion of a driven 1-dimensional system of N oscillators in the presence of 
damping has been well analyzed in the framework of Newtonian mechanics. 
However, it has not had the same attention in higher dimensions. The basic 
problem lies in the fact that since the motion of an oscillator in a lattice 
is influenced by its nearest neighbors, the equations of motion are coupled 
amongst the perpendicular directions, viz., x, y and z components. However 
under the approximation that the perturbations are small, the motion can 
be decoupled in the individual directions. I discuss the problem of free 
oscillations in the first section and introduce damping and driving terms in 
the following section. 
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2 Undamped Natural Oscillations in 2D 



In this section, I will discuss the effect of a perturbation on a 2-dimensional 
lattice of bodies identical in mass. The lattice under consideration is a finite 
(NxM) structure. I deal with systems under natural oscillations, viz., under 
the absence of damping and external forces. Every oscillator is connected 
with its closest neighbor with a linear spring. For simplicity, the lattice is 
understood to lie in the first quadrant of the regular cartesian co-ordinate 
system. For integral values of i and j, the ordered pair denotes a body 
belonging to the i th column and j th row. The initial distance between any 
two closest neighbors is equal to the unstretched length of the spring, a. So 
the position vector fij of an oscillator can be written as 

fij = %■ ax + j ■ ay 

On a small perturbation, the new position vector of an oscillator r^j 
is given by 



1,3 



i-ax + j -ay + r'ij (1) 



The equation of motion for a body not on any of the boundaries can be 
written as 

cPn ■ 

m ~^f = F (i-i,MhJ) + F (i+i,j),(i,j) + F (i,j-i),(i,j) + F (i,j+i),(i,j) ( 2 ) 

where uiq 2 is the spring constant and m is the mass of the oscillator. Note 
that the convention used here is that of Goldstein; F^ v represents the force 
on the particle v by the particle fx. can thus be expressed as 

F (i-i,j),(i,j) = -^o 2 [ I fi,j ~ rt-ij I -a] y _ (3) 
For small perturbations 

r i,j ~ r i-l,j r i,j ~ r i~l,j 

r^j 

I r i,j ~ r i-l,j I a 

Using (1), 



yj {xi,j ~ Xi~i,j) 2 + (yij - yi~i,j) 2 

^+^- I S-ij) 2 + ^-ii J ) ! 
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For small perturbations, 0(x /a ) 



J 1 + Ht_^=a 



On binomial expansion of the square-root term under the above approxima- 
tion and using (3), force can be written as 



-» r 7 x i— 1 7 "i 

fi j - Ti-t j -a = a 1 -I - - a 

a a J 



Expressing the position vector in terms of its components and using (1), 

a a a a a 

Neglecting terms of order 0(x' 2 /a 2 ), 0(y' 2 /a 2 ) and 0(x'y' /a 2 ), (4) can be 
expressed as 

The other terms in (2) can be evaluated similarly and the equation of motion 
for can be written as 

d 2 V' ' 

m —[p - = - w o[( 2 4j - x 'i-i,j - x 'i+i,j) £ + ( 2 y'i,j - y'i,j-i - y'i,j+i)y] 

One can see that the equation of motion is decoupled into two individual 
motions, about x and y such that, 



m ~d^~ = = [ 4 > j ~ Xi ~ lj ~ X ' l+1 ' j] ® 

d 2 Ui,j _ d 2 y i j 2 r / / / i / c \ 



Once the motion is decoupled in the individual directions, the problem re- 
duces to the familiar system of N coupled oscillators in one dimension pQ. 
Assuming a test solution of the form, 

x\ j = Aij cos(u)t) (7) 
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and plugging it in the equation of motion (5) one gets 

A i+1J + Ai- ltj -w 2 + 2 W 2 

Since the ratio of the amplitudes is constant, one can assume Aij to be of 
the form, 

A i:j = Csm(i6) cos(j(f)) (9) 
On substituting the above expression in (8), 

w = 2w sin(^) (10) 



The values of 9 and <p can be found by imposing the boundary condition 
that the ends of the lattice are held fixed, hence 

Xoj =0 Sz XM-1J = 
x ifi = & x it N-i = 

On considering the boundary values, it is seen that 9 and (f> need to obey 
sin ((M - 1)9) =0 & sin ((N - 1)9) = 

Hence, 



mir 
M - 1 



where m = 0,1, ... ,M - 1 (11) 



Ti 7T 

(f>=jj — j- tuyere n = 0,1,. .. ,7V - 1 (12) 

The indices, m and n correspond to the particular normal mode of the 
system. One can see that for a N x M lattice of oscillators, there exists 
MN normal modes. Thus the motion of the system can be described by(up 
to an arbitrary phase factor) 



M-1AT-1 , v , , 

i=0 j=0 ^ ' ^ ' 



to carries an index m, since by (11) the value of 9 for different values of m 
are different. So by (10) the frequency is seen to depend on m. One can 
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similarly treat the motion in y with the boundary conditions that the ends 
are held fixed and obtain(up to an arbitrary phase factor), 



uj n = 2ujq sin - (14) 



M-lN-l , v , s 

y' = E E Ann sin (i-^^ J sin (i^^y ) cos(w n t) (15) 

i=0 j=0 ^ ' ^ ' 

In the above expression 6> and (f> take the same values as (11) and (12). The 
values of C mn and D mn can be found by imposing initial conditions. 

2.1 Undamped Natural oscillations in 3D 

Natural oscillations in 3-Dimensions can be treated in a similar way. For a 
M x N x P lattice under the boundary conditions that the ends are held 
fixed, following the convention developed in the parent section(up to an 
arbitrary phase factor), 



m — i iv — i r — i x s ✓ v x v 

= XI EE C ranpSin (*^37 ) sin UjfZTl ) ^ l^F^T ) COS ( W "^) 
i=0 j=0 fc=0 ^ ' ^ ' ^ ' 

M-iJV-iP-i f m-K \ ( nir \ ( it \ 

y = E E E sin ( * ) sin ( J ' iv~T ) sin ( ^P^T ) cos ^*) 



j=0 j=0 fc=0 
M-l AT-1 P-l 



/ \ - \ - v - / mir \ . / . nir \ . I pir \ , . 



i=0 j=0 k=0 

The frequencies are given as 



where 



9 <f> ip 

u m = 2w sm(-) w n = 2w sin(-) cj p = 2w sm(— ) (16) 



7717T 

where m = 0, 1, . . . ,M - 1 (17) 



M - 1 
n7r 

N — 1 



where n = 0, 1, . . . , JV - 1 (18) 



^ = -21- where p = 0,1,... ,P -1 (19) 
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3 Forced oscillations with Damping 



In this section, I aim to study the motion of a system in the presence of a 
first order velocity dependent damping force and an external periodic force. 
I shall once again consider a M x N lattice as discussed in the previous 
section. If one tries to apply the approach similar to the preceding section, 
it is seen that the ratio of the amplitudes (8) take up time dependent values, 
hence one cannot assume the amplitudes to be of the form (9). If the system 
is under an external driving force given as 



h3 



(AijX + Bijy)e 



-iuiQt 



(20) 



where Aij and Bij are the amplitudes of the external force. The damping 
force is taken to be a velocity dependent force to first order, given as 



ampij 



^ dt X + Ji ' j dt V 



(21) 



dij and fcj are the damping coefficients. One can revisit (1) and (3) to 
incorporate the given damping and driving force so as to write the equations 
of motion of the oscillator as 



m 



m 



dt 2 

d 2 yij 



dt 

d Vi, 



+ u (2x' i j - — x i+i,j) — Aije 



-iu>ot 



dt 2 ' JhJ dt ' ~»^3 

For a particular value of j = jo, (22) can be represented as 
K X + D-X + M • X = Fe-™ * 



(22) 
(23) 

(24) 



Where K, M, D & F are the mass, dissipation, spring constant and driving 
force matrices and uq is the frequency of the driving force. The matrix 
representations are given as 



X = 



x i,j 

V x M-l,j J 



D 



/ do,j 



v '■■ 




dij 
d 2 j 



\ 



6 



/ A o,j \ 



2,7 



K 



U 2 Q 
-UJ 2 





-U) 2 

2w 2 
-w 2 






-u 2 o 



V 



/ 



For the sake of simplicity, as mentioned before all the oscillators are assumed 
to be of the same mass, m hence the mass matrix is an M x N diagonal 
matrix with the masses of the oscillators as its diagonal entries. The solution 
of the above equation is of the form [2] 

X = (K - iuj Q T) - oj^M)- 1 ■ Fe~ iuJot + X /ree (25) 

where Xj ree is the solution for homogeneous conditions imposed on (24), 
viz., under the absence of external driving force. 



X 



free 



-iuj n t 



(26) 



Hence the problem practically reduces to finding the coefficients C n in terms 
of the initial boundary conditions and solving the eigenvalue problem for the 
operator. The coefficients can be expressed in terms of the initial conditions 
as 

X„ • F 



x = ^x„ 



-iujot 



X n ■ [(u n M + %D) ■ X^ ree (0) + iM ■ X^ ree (0)]( 



— iuJnt 



(27) 



Hence the motion along the x axis is completely determined once the initial 
conditions are known. Similar treatment of the equation of motion in the y 
direction, (23) gives the solution for the motion along the y axis. One can 
also extend this approach to 3-dimensions. 



4 Conclusions 

It is seen that once the motion is decoupled in the individual directions the 
motion of the system can be expressed as an expansion amongst the Normal 
modes. I aim to study the non linear dynamics of the system |3j [4 J |5j,[6j to 
solve for the motion in the presence of non- linear terms, viz., 0(x' 2 /a 2 ) and 
higher. 
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